Convergence of adaptive mixed finite element method for convection-diffusion-reaction equations

We prove the convergence of an adaptive mixed finite element method (AMFEM) for (nonsymmetric) convection-diffusion-reaction equations. The convergence result holds from the cases where convection or reaction is not present to convection-or reaction-dominated problems. A novel technique of analysis is developed without any quasi orthogonality for stress and displacement variables, and without marking the oscillation dependent on discrete solutions and data. We show that AMFEM is a contraction of the error of the stress and displacement variables plus some quantity. Numerical experiments confirm the theoretical results.Comment: arXiv admin note: text overlap with arXiv:1312.645

Second Order Reliability Method for Time-Dependent Reliability Analysis Using Sequential Efficient Global Optimization

Reliability depends on time if the associated limit-state function includes time. A time-dependent reliability problem can be converted into a time-independent reliability problem by using the extreme value of the limit-state function. Then the first order reliability method can be used but it may produce a large error since the extreme limit-state function is usually highly nonlinear. This study proposes a new reliability method so that the second order reliability method can be applied to time-dependent reliability analysis for higher accuracy while maintaining high efficiency. The method employs sequential efficient global optimization to transform the time-dependent reliability analysis into the time-independent problem. The Hessian approximation and envelope theorem are used to obtain the second order information of the extreme limit-state function. Then the second order saddlepoint approximation is use to evaluate the reliability. The accuracy and efficiency of the proposed method are verified through numerical examples

A Bayesian Approach to Recovering Missing Component Dependence for System Reliability Prediction via Synergy Between Physics and Data

Predicting system reliability is often a core task in systems design. System reliability depends on component reliability and dependence of components. Component reliability can be predicted with a physics-based approach if the associated physical models are available. If the models do not exist, component reliability may be estimated from data. When both types of components coexist, their dependence is often unknown, and the component states are therefore assumed independent by the traditional method, which can result in a large error. This work proposes a new system reliability method to recover the missing component dependence, thereby leading to a more accurate estimate of the joint probability density (PDF) of all the component states. The method works for series systems whose load is shared by its components that may fail due to excessive loading. For components without physical models available, the load data are recorded upon failure, and equivalent physical models are created; the model parameters are estimated by the proposed Bayesian approach. Then models of all component states become available, and the dependence of component states, as well as their joint PDF, can be estimated. Four examples are used to evaluate the proposed method, and the results indicate that the proposed method can produce more accurate predictions of system reliability than the traditional method that assumes independent component states

Active learning with generalized sliced inverse regression for high-dimensional reliability analysis

It is computationally expensive to predict reliability using physical models at the design stage if many random input variables exist. This work introduces a dimension reduction technique based on generalized sliced inverse regression (GSIR) to mitigate the curse of dimensionality. The proposed high dimensional reliability method enables active learning to integrate GSIR, Gaussian Process (GP) modeling, and Importance Sampling (IS), resulting in an accurate reliability prediction at a reduced computational cost. The new method consists of three core steps, 1) identification of the importance sampling region, 2) dimension reduction by GSIR to produce a sufficient predictor, and 3) construction of a GP model for the true response with respect to the sufficient predictor in the reduced-dimension space. High accuracy and efficiency are achieved with active learning that is iteratively executed with the above three steps by adding new training points one by one in the region with a high chance of failure

High-Dimensional Reliability Method Accounting for Important and Unimportant Input Variables

Reliability analysis is a core element in engineering design and can be performed with physical models (limit-state functions). Reliability analysis becomes computationally expensive when the dimensionality of input random variables is high. This work develops a high-dimensional reliability analysis method through a new dimension reduction strategy so that the contributions of unimportant input variables are also accommodated after dimension reduction. Dimension reduction is performed with the first iteration of the first-order reliability method (FORM), which identifies important and unimportant input variables. Then a higher order reliability analysis is performed in the reduced space of only important input variables. The reliability obtained in the reduced space is then integrated with the contributions of unimportant input variables, resulting in the final reliability prediction that accounts for both types of input variables. Consequently, the new reliability method is more accurate than the traditional method which fixes unimportant input variables at their means. The accuracy is demonstrated by three examples

Uncertainty Quantification by Convolutional Neural Network Gaussian Process Regression with Image and Numerical Data

Uncertainty Quantification (UQ) plays a critical role in engineering analysis and design. Regression is commonly employed to construct surrogate models to replace expensive simulation models for UQ. Classical regression methods suffer from the curse of dimensionality, especially when image data and numerical data coexist, which makes UQ computationally unaffordable. In this work, we propose a Convolutional Neural Network (CNN) based framework, which accommodates both image and numerical data. We first transform numerical data into images and then combine them with existing image data. The combined images are fed to CNN for regression. To obtain the model uncertainty, we integrate CNN with Gaussian Process (GP), which results in the mixed network CNN-GP. The simulation results show that CNN-GP can build accurate surrogate models for UQ with mixed data and that CNN-GP can also provide the uncertainty associated with the model prediction